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Abstract 

Supersymmetric quantum mechanics has many applications, and typically 
uses a raising and lowering operator formalism. For one dimensional problems, 
we show how such raising and lowering operators may be generalized to include 
an arbitrary function. As a result, the usual Rodrigues’ formula of the theory 
of orthogonal polynomials may be recovered in special cases, and it may other¬ 
wise be generalized to incorporate an arbitrary function. We provide example 
generalized operators for several important classical orthogonal polynomials, 
including Chebyshev, Gegenbauer, and other polynomials. In particular, as 
concerns Legendre polynomials and associated Legendre functions, we supple¬ 
ment and generalize results of Bazeia and Das. 
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Introduction 


The concepts of supersymmetry from quantum field theory have led to diverse ap¬ 
plications in many other areas of physics including statistical, condensed matter, and 
nuclear physics, and quantum mechanics [Hiiniiig. The supersymmetric scheme per¬ 
mits to combine fermionic and bosonic degrees of freedom. Supersymmetry presents 
a hierarchy of Hamiltonians, and for nonrelativistic quantum mechanics provides fac¬ 
torizations of the Schrodinger differential equation. Factorization of the Hamiltonian 
is done with operators incorporating the superpotential W, and W satishes Ricatti 
equations with the partner potentials being the inhomogeneous term. 

The factorization of Hamiltonians is typically developed via raising and lowering 
operators, each being of hrst order. In this paper we focus on the one dimensional 
Schrodinger equation and describe how an arbitrary function, denoted h(x), may be 
incorporated into these operators. This procedure allows known Rodrigues’ formulae 
to be recovered in short order for even the simplest choices of h(x), and otherwise 
gives a way to develop generalized Rodrigues’ formulae and other identities. In this 
Introduction we provide illustrations for the quantum simple harmonic oscillator, and 
the associated Hermite and Laguerre polynomials. We then concentrate on the Leg¬ 
endre differential equation, thereby generalizing many of the results of [1]. We also 
supplement |1] by giving proofs of certain identities. Furthermore, we provide a gen¬ 
eralization to Gegenbauer polynomials with parameter A. Our method is applicable 
in a great many contexts, and we briefly mention other examples at the end. 

We let L^{x) denote the nth Laguerre polynomial, Hn{x) the nth Hermite poly¬ 
nomial, (z)n = z{z + l) ■ ■■{z + n — 1) =T{z + n)/T{z) the Pochhammer symbol, with 
F the Gamma function (e.g., mmm)- We recall that Hermite polynomials are 
special cases of Laguerre polynomials, such that 

H 2 nix) = i-ir2^^n\L-^/^{x^) and H 2 n+iix) = {-ir2^^+^n\xLl/^x‘^) (1.1) 
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(e.g., HU, p. 1037). Both of these types of polynomials are instances of the continent 
hypergeometric fnnction iFi. 

We work with dimensionless qnantities thronghont, so that for the qnantnm har¬ 
monic oscillator the raising and lowering operators may be taken as For 

the raising operator, we make the ansatz 


, d „ d . 

a — X ~ ^ + h, 

dx dx 


wherein the fnnctions /i and g 2 are to be determined in terms of the arbitrary fnnction 
h. Letting act on a differentiable fnnction /, and then eqnating coefficients of / 
and /', we hnd that, aside from mnltiplicative constants. 


g 2 {x) = exp 




and 

(x — h)dx . 

It goes very similarly for generalizing a~, so that with a sign change the roles of fi 
and g 2 are reversed. The original operators do not depend npon the energy level n 
and accordingly neither do the generalized operators with h in general different from 
X. The well known eqnal energy spacings follow. 

As usnal, the gronnd state obtains from aipo = 0 and 'ipo{x) oc exp(—x^/2). Acting 
snccessively on the gronnd state wavefnnction with the raising operator provides the 
Hermite fnnctions. If h = 0, then a'^ = —By mnltiplying by 
and imposing a sign convention, even in the simplest case of taking h = 0, the nsnal 
Rodrignes’ formnia for Hermite polynomials is recovered. 

For the Lagnerre differential eqnation 


fii^) = 


92{x) 


= — exp 


xu" + {a — X + l)u' -|- nn = 0, 


( 1 . 2 ) 
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raising and lowering operators may be taken as 

. d A d 

A+ = X- - X + a + n, = -x -—h n, 

ax ax 

so that the differential equation assumes the forms 

A+A_L'^{x) = n{n + a)L"(a;) 


(1.3) 


and 

A_A+Ll_^{x) = n{n + a)Ll_^{x). 

Here we have applied recurrence relations m, 8.971.3) so that = nL°{x) 

and A_L'^{x) = (n+a)i^“_i(x). If we put A^ = fi-^g 2 +h, then we find that, omitting 
a multiplicative constant, 


g2{x) = 


X^+rig ^gxp 

1- 

1 


J X 


and 

f^{x) = = x^-^-'^e^exp [ [ ^dx . 

g2{x) ' [j X \ 

The correspondence to the usual Rodrigues’ formula T"(x) = (e“’^x”+“) 

is then evident again for h = 0. The Laguerre polynomials appear, for instance, in 
the wave function for the hydrogen atom. 


Generalized raising and lowering operators for the Legendre equations 


We now focus on the Legendre differential equation (x = cos^, —1 < x < 1) 

Pn{x) = 0 , ( 2 . 1 ) 


d , n , d 


and the associated Legendre differential equation 


d , n , d 
— (x^ - 1):^ + 


dx 


dx 1 — x^ 


— n{n + 1) 


Pn{^) = 0, 


( 2 . 2 ) 
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with \m\ < n. The associated Legendre functions are related to the Legendre poly¬ 


nomials via 


Jra pm(^\ 

P^(X) = (1 - W 

^ ^ ^ ’ dx^ 


Note that for m odd, P^{x) is algebraic, not polynomial, in x. 
Taking the raising operator for (2.1) to be 


d 


d 


Rn = (x - 1)— + nx = fi^g 2 + h, 
dx dx 

we hnd, multiplicative constants henceforth being omitted, 

92{x) = (x^ - l)”/^exp 


-1 

I 


L J 

' x^-l 


and 


— 1 


/i(a^) = ^T = (^ -1) exp 

92[X) 

As a brief example, when h{x) = ax^ + hx + c, so that 


h{x) 
x"^ — 1 


-dx 


ax^ + hx + c. 1., . \ 1 \ \ 

-r- dx = ax —(a -|- 0 -|- c) lull — x) -|- (o — a — c) lull -|- x), 

x^ — I 2 


the exponential factor becomes 

h{x 


exp 


— I 


dx 


= e“^(l - . 


Proceeding similarly with the lowering operator for (2.1), 


Ln = -I) +nx = - 92 -^ fi + h, 


we determine that 


Mx) = (x^ - l)'-"''2exp 


h{x) 
x"^ — I' 


-dx 


and 


— I 


92 ix) = = (x - Ip exp 


-1 

I 

f w ' 

/ - dx 

L J 

' a;2-l 


(2.3) 
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Generalizing the m-raising operator for the associated Legendre eqnation (2.2), 

we have fi{x) = (1 — /g 2 {x) and 


d mx d 

+h, 


g^(x) = (1 -a;2)-"*/2exp 
Similarly for the m-lowering operator 


h{x) 


(1 _3.2)1/2 


dx 


L =- — (i-x^Y^^ + _ — _ 

dx^ ’ ^(l-a;2)V2 


— — 92 -^fi + h, 
dx 


we determine that 


fYx) = (l-a;2)("*+^)/'exp 


h{x) 


(1 _a.2)l/2 


dx 


with g 2 {x) = (1 — 

These operators satisfy i?„Pn_i(a;) = nPn{x), Ln+2Pn{x) = nPn-i{x), RmPn{x) = 
and LmPn^^{x) = [n(n +1) — m(m +l)]P^(a;) = (n —m)(m + n +l)P™(a:). 
For h = 0 these resnlts reduce to the raising and lowering operators of (10) and (27) 

of [1]. 

From the generalized raising operator (2.3) we find that 


(^^ — 1) 

Pn{x) = --exp 


_l(a;2_i)n/2exp 

dx 


[ h{x) 

x"^ — 1 

+h{x)- - -{x^ - 1) 

^ ^ dx'^-^ ^ ^ 



f h{x) 

d^- 

i J 

' x'^ — 1 

dx^ 


2 1\n-l 


n -1 


— (X^-I) 

7—1 ^ ^ 


(2.4) 


This leads to the equation 2"' ^n\Pn{x) = Rn'P' ^{n — l)!P„_i(a:). The relation (2.4) 
is a generalization of (16) of |1] for nonzero functions h{x). 

We let 

h{x) 


e[x) = exp 


x^ — 1 


(2.5) 


Then we may arrive at families of generalized Rodrigues’ formulae in the form 


n!P„(a:) = (a;2-l)W2[e(^)]n-i^ 

dx 


(^2 _ 1)3/2^ 

dx 




{x ^ - 1)V^ 
e{x) 


+ F[ft|. 
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wherein F, being 0 for h = 0, depends upon powers of h and its derivatives with re¬ 
spect to X. The term of F with the highest order derivative of h is — [x) 
and the term of the highest power of h is (—l)"'(n — l)\xh'^~^{x). 


Generalization to Gegenbauer polynomials 

One of the generalizations of Legendre polynomials is the Gegenbauer (or ultra- 
spherical) polynomials in turn a special case of the Jacobi polynomials P^'^\ 


CGx) = 


(2A), 




(A -l-1/2), 

These polynomials satisfy the ordinary differential equation 


— l)y" + (2A -I- l)xy' — n(2A -|- n)y = 0, 


and the case A = 1/2 recovers the Legendre polynomials. A pair of lowering and 
raising operators in n for the Gegenbauer polynomials is given by 

(1 - + nx C^{x) = {n + 2X- l)G^_i(x), 


Putting 


we hnd that 


and 


(1 - - (n - 1 + 2A)x C^_^{x) = -nC^{x). 


G+ = (1 - x^)^ - (n - 1 2A)x = /i^ 5'2 + K 


dx 


= n - .,,2yn-l+2A)/2 


dx" 


92ix) = (1 - X ) 


exp 


h{x) 
x^ — 1' 


-dx 


-1 

1 

f J 


T—1 

1 


f^[x) = = (1 -x2)(3-G/2-Agxp 

g2[x) 

With Gq (x) = 1, via iteration of G^(x) = — .;^(x) for h = 0 we determine that 

n—1 


C^^X) = ^^(1 - a;2)(3-n)/2-A^ 
TZ. CtOC 


( 1 _ ^ 2 ) 3/2 

dx 


{1-xT- 
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For comparison, the usual Rodrigues’ formula for Gegenbauer polynomials is 


C^(r) = ( 2)"'r(n + A)r(n + 2A) _ ^ 2 '>i/ 2 -a_^/i _ 2 '>n+A-i /2 

^ n\ r(A)r(2n + 2A) ^ ^ ’ 

From the latter formula we determine that 




( 2 )n-i Y(n + X- l)F(n + 2A - 1) 


n\ 


! F(A)F(2n + 2A-2) 


X(1 _3.2)(3-n)/2-A_^^^ _^2^n/2_^(^ _^2)n+A-3/2_ 

(] r 1 — -L 


dx^' 


A Legendre polynomial identity 


The following identity, useful for deriving a Rodrigues’ type formula, was given 
without proof in |1]. We present three proofs of this identity, the hrst using corre¬ 
spondences with the Legendre polynomials, and the other two employing the product 
(Leibniz) rule. For n>2, 

Jn J7i 2 

For the hrst proof, we begin by making various identihcations of the left side of 
the identity (3.1) in terms of Legendre polynomials Pn{x) and associated Legendre 
functions P^{x) Plain]. From the usual Rodrigues’ formula 

1 

the left side of (3.1) is 

2”-‘(n - l)!(a:= - l)^P„_,(x) = -2”-‘(n - 1)!(1 - xY'"P^_i(x), 

where we have used Pl{x) = (1 — x‘^Y^‘^-£:Pe{x). From recurrence 12.5.8(f) of [5] we 
also have 

- l)^Pn-l{x) = 2'^~^n\[Pn{x) - xPn-l{x)]. 
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2"-^(n- l)!(a;2 


(3.3) 









The right side of (3.1) is the integral of the polynomial 


n-l 


d 

- ^)n\Pn-i{x) 

By integrating the recurrence 12.5.8(b) of [5] by parts we obtain 


(n + 1) / Pn{x)dx = -Pn-l{x) + xPn{x), 


(3.4) 


and 


(n + 1) / Pn{x)dx = -Pn-l{x) + xPn{x) + P„_i(0). 
Jo 


The right side of (3.1) is, using (3.4), 


{n - l)n j {x - I)*" dx = 2^ (n - l)n! J Pn-i{x)d. 

= (n - l)!(n - 1 ) 2 ”“^[-P„_ 2 (a;) + xPn-i{x)] 

= n\2'^~^[Pn{x) - xPn-i{x)]. 


X 


(3.5) 


In the last step we have used recurrence 12.5.8(a) of j3] in order to eliminate the 
Pn- 2 {x) term. The equality of (3.3) and (3.5) verifies identity (3.1). □ 

Remark. The expressions (3.3) or (3.5) verify that (3.1) evaluates to 0 when 
X = ±1. 

A second proof of (3.1) may be based upon the product rule and the factorization 
— 1 = (x — l)(a: + 1). We then obtain for the left side of (3.1) 

(" ) (1 - ")'(1 - ”)»-<(^ - l )'(“^ + 1 )”^'- ( 3 - 6 ) 


while the right side is given by 


£=0 


n—2 


n[n 


i)(-i)”E 


i=0 

n—1 


= n(n- l)(-l)”^ 


e=i 


n — 2 


n — 2 


(1 - n)n_^_ 2 (l - n)i{x - lY"^^{x + 1) 




(1 - n)n-i-i{l - n)i_i{x - lY{x + 1)'" (3.7) 


9 





Since in (3.6), (1 — n)n = 0 unless n = 0, the £ = 0 and i = n terms evaluate to 0. 
Comparing (3.6) and (3.7) we need to verify the equality 


- 1) (1 - n)n-£-i{l - n)i_i = Q {1 - n)e{l - n)n-e- 

With the ratios (1 — n)n-e/{l — n)n-e-i = —i and (1 — n)£/{l — n)i_i = i — n, this 
equality reduces to the readily conhrmed statement that 


n{n — 1) 


n — 2 
- 1 


nl 


[n - t t = 


{i - i)\{n - i - ly: 


□ 

For a third proof, we hrst binomially expand each side of (3.1), and then differ¬ 
entiate term by term: 


^ ) /rt — 


dx 


dx 


e=o 


\n-i-lJ2i 


and 


n—1 




£=0 


n — 1 




2l-n+2 


d” ^ /n — 1 


dx"" 


k=0 ^ 


n—1 


-2(-irY, 


k=0 


n — 1 
k 


{-l)-’^-^k{l-2k)n-iX 


2k—n 


Therefore, the left side of (3.1) is given by 


n—1 


k=0 


k 




2k—n+2 _ ^2k—n\ 


= -2(-i)^ 


n—1 

E 

. fc =0 


n — 1 




2k—n+2 


n—2 


E (l + !) + 1)(-1 - 2k)„_,x 


2k—n-\-2 


/c=0 
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Thus, equating the coefficients of (3.1) is valid if the following equality is 

verihed: 





= n{n — 1) 


n — 1 


(-l)”-‘-V(l - 2k)n-3. 


(3.8) 


We have the ratios 


(1 - 2k)n-3 _ 1 

{-l-2k)n-i 2k{2k + l)' 


(1 - 2k)n-i 
(1 — 2k)n-3 


{2k — n + l){2k — n + 2), 


and 

(1 — 2A;)„_i {2k — n + l){2k — n + 2) (fc+i) n — k — 1 

{-l-2k)n-i ^ 2k{2k + l) ’ ^ k + 1 

Then indeed (3.8) reduces to the equality 


(2A:-n + l)(2A: 

WTT 


n + 2) 


+ 2(n 


k-1) 


n{n — 1) 
2A: + 1 ■ 


□ 


Remark. It may be shown that (3.1) is equivalent to the three-term identity 

Jn—1 

(^^ - = (^ + l)(^ - - l)"■^ 

which in terms of Legendre polynomials is the equality 

2r 

(X^ - l)Pn-l{x) + - - -[Pn-2ix) - Pnix)] 

2n — 1 


{n -I- l)(n — 2) 
1 - 2n 


2n — 3 


[Pn-lix) - Pn-3ix)] + 


2n -I- 1 


[Pn-lix) - Pn+lix)] 


Discussion 


For the Chebyshev polynomials of the second kind Unix) natts], a set of lowering 
and raising operators is given by 


a 2\ ^ 

— X ) -—\-rix 
ax 


Unix) 


in + l)Un-iix), 
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and 


(x^ - 1)— + (n + 2)x Unix) = (n + l)Un+i(x). 


With [/+ = (x^ — 1)^ + (n + 2)x = f 1-^92 + h, we may replace n —)■ n + 2 in the 


result of (2.3), giving 


— ( 1 \n/2+l _ 


g^[x) = ix - 1 ) 


eix) ’ 


hix) = ^ I I = - 1 ) '"^‘^eix), 

92ix) 


with e{x) as given in (2.5). Taking the special case of h = 0, with Uq = 1 and the 
iteration nlUnix) = U^_-^U ^_2 ■ ■ ■ U^UqUq, we obtain 

7 r g 1 

nlUnix) = (x^ — l)d-»^)/2 (^ 2;2 _ l'^3/2 ^^2 _ 

dx I dx\ 


This expression may be contrasted with the trigonometric dehnition 


UJx) = 


sin[(n + 1) cos ^ x 
sin(cos“^ x) 


and the usual Rodrigues’ formula 

JJ ^ (-l)^(n + l) dU _ ^2^|n+l/2 

^ yr^(2n + l)!!da;-^ ^ 

where (2n + 1)!! = (2n + l)(2n — 1) ■ ■ - 5 ■ 3 ■ 1. Using this formula and Unix) = 
ix) gives the expression 


The Chebyshev polynomials of the hrst kind Tnix) satisfy the differential equation 
(1 — x‘^)u'' — xu' + ri^u = 0 and may be dehned from Tnix) = cos[ncos“^ x\. Raising 
and lowering operators such that T^T^ix) = Tm±iix), are given by, for m > 0, 


= =F 

m. * 


(1 — x^) d 


m dx 


- — X. 
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Putting T+ = /i^fi '2 + h, we find that 


g^(x) = (1 -x2)™/2exp 


m 


h{x) 

1 — x‘^ 


:dx 


and 


h{^) = = -1(1 - a;2)—/2+1 

mg 2 {x) m 


exp 


—m 


h{x) 

1 — x‘^ 


-dx 


We have Tq = 1, Ti{x) = x, and Tm{x) = T+_^T +_2 • • ■Tj+Ti(x), giving for h = 0 

-| n—2 


TJx) = 


LlUfl _ 3.2')(3-n)/2_^ 
{n-l)r ^ dx 


(1 - x2)3/2 


d 

dx 


x{i — x'^y^'^ 


In comparison, the usual Rodrigues’ formula for Tn{x) is 

{-lyyi-x^ d 


Tn{x) = 


(l_^2)n-l/2_ 


(2n — 1)!! dx^‘ 

Using this formula and Tn{x) = Ty_^Tn-i{x) gives the expression 


Tn{x) = 


■1)” _ 3.2)(3-n)/2_^^^ _ _ 3.2)n-3/2_ 

rirjr> ri rv^Tl 1 


(2n — 3)!!(n — 1) ^ ^ dx'^ ^ dx'"' 

With a change of variable x = cosy, the raising and lowering operators for and 
Un have been identified as elements of a sm(1, 1) algebra for the infinitely deep square 
well potential. They were also used to study the temporally stable coherent states 
for the infinite well and Poschl-Teller potentials [3]. 

For the radial Coulomb problem, a raising operator may take the form 

with angular momentum quantum number i (cf. [El [HI E]). Then we hnd fi{r) = 
l/g 2 {r) and 


g 2 {r) = exp 


h{r)dr 


Similarly, for a three dimensional harmonic oscillator with raising operator 

, d r £ + 1 j, , . d , . ,,, 
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we have /i(r) = llg 2 {r) and 


g 2 {r) = exp — j h{r)dr 


For these two (scaled) radial problems, the wave function contains products of powers 
of r, a decreasing exponential factor exp(—r/2), and Laguerre polynomials. 

Many more applications of our approach are possible. In some more generality, 
suppose, for example, that a lowering operator is given as 

L = + h, 


wherein a and b are specihed functions. Then we determine that 


fi{x) = a{x)exp 


^ ] \ f ^ ' 

ax exp / ax , 
a[x) \ \_J a[x) \ 


a{x) 

92{x) = -7^ = exp 
fi{x) 


^ ] \ f ^ 

ax exp — / ax 
aix) J aix) 


Indeed, we have produced generalized lowering and raising operators and a general¬ 
ized Rodrigues’ formula for arbitrary Jacobi polynomials Pji^’^\x), and this will be 
reported elsewhere. 

Another direction for obtaining Rodrigues-type formulae proceeds as follows, il¬ 
lustrated for the Laguerre polynomials of the Introduction. Making the change of 
variable x = P, the differential equation (1.2) becomes 

+ + ^ + 4:ny{r^) = 0 , 

ar^ r \ 2 J dr 

and the raising and lowering operators (1.3) transform to 

. r d , , r d 

The generalized raising operator A+ = fi{r)^g 2 {r) + h{r) has 


g 2 {r) = exp -2 


dr 
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and fi{r) = r/[2g2{r)]. By iteration with h = 0 and Lg(r^) = 1 we obtain 

L«(r2) = ]_^er^^-2in+a)d / ^3 d V ' _ 

n\ 2” dr \ dr J 

Such relations may be all the more important because of the connection (1.1) to 
Hermite polynomials when a = ±1/2, so that 

H,4r) = (-l)”2'V-“e"± 

and 

In conclusion, we have generalized raising and lowering operators that are so use¬ 
ful for quantum mechanics. These then lead to generalized Rodrigues’ formulae and 
other identities. The generalized operators have been illustrated for the Chebyshev, 
Gegenbauer, Hermite, Laguerre, and Legendre polynomials. Especially the latter 
three polynomials appear in important quantum mechanical and other problems, in¬ 
cluding the harmonic oscillator and Coulomb problem. Moreover, when combined 
with Gaussian factors, Hermite and Laguerre polynomials are very useful for describ¬ 
ing optical helds. For a variety of problems with spherical symmetry, the Legendre 
polynomials and associated Legendre functions appear. We have provided discussion 
specihc to the latter functions that generalize previous results. Given the existence 
of generalized Rodrigues’ formulae, there is also the possibility to establish general¬ 
ized Schlaefli integrals. The latter contour integrals also play a role in the theory of 
orthogonal polynomials. 
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